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Abstract. We study a magnitude-limited sample of 10 gamma-ray burst (GRB) host galaxies with known spectroscopic 
redshifts (0.43 < z < 2.04). From an analysis of the spectral energy distributions (SEDs), based on published broad-band 
optical and near-infrared photometry, we derive photometric redshifts, galaxy types, ages of the dominant stellar populations, 
internal extinctions, and ultraviolet (UV) star-formation rates (SFRs) of the host galaxies. The photometric redshifts are quite 
accurate despite the heterogeneous nature of the sample: The r.m.s. errors are cr(z) = 0.21 and cr(Az/(l + Zspec)) = 0.16 
with no significant systematic offsets. All the host galaxies have SEDs similar to young starburst galaxies with moderate to low 
extinction. A comparison of specific SFRs with those of high-redshift galaxies in the Hubble Deep Fields shows that GRB hosts 
are most likely similar to the field galaxies with the largest specific SFRs. On the other hand, GRB hosts are not significantly 
younger than starburst field galaxies at similar redshifts, but are found to be younger than a sample of all types of field galaxies. 
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1. Introduction 

The association of Gamma-Ray Bursts (GRBs) with collapsing 
massive stars has been inferred through many observations dur- 
ing the past 7 years, notably through supernova bumps in the 
afterglow light curves or localisation s of the afterglow close 
to star forming regions ( Castro-Tirado & Gorosabel 1999; 
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Bloom et al. 1999; Fruchter et al. 


1999a; 


Galama et al. 


1998; 


Kulkarni et alJ Il998j 




Holland & Hiorth 


1 19991 iFvnbo et alJ 



genitors, given their small positional offsets relative to their 
hosts, are likely massive stars. Evidence that at least some 
long-duration GRBs are associated with supernovae came with 
the observations of the afterglow of GRB 030329, which 
showed distinct spe ctral supernova features JStanek et alJ2003l 
iHiorth et alJ2003bh . 

Underlying host galaxies have been found in all cases 
where GRBs were localised with sub-arcsecond precision. 
Currently, the sample consists of ~35 s uch hosts with known 
redshifts in the range 0.1055 < z < 4.5 JProchaska et all2 004: 
lAndersen et alJ200ol). The reported magnitudes of the hosts are 
20.4 < R < 30( lMalesanietalJ2004[|JaunsenetalJ2003l) . 

The life time of the massive stars believed to produce long- 
duration GRBs is of the order of a few Myr. If the host galaxies 
are indeed forming such massive stars this should be reflected 
in their spectral energy distributions (SEDs) which, in the ab- 
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sence of reddening, should reveal a blue continuum. Moreover, 
their star-formation rates (SFRs) should be large. The inte- 
grated SFR of a galaxy can be found from measurements of 
Ha or [Oil] line fluxes, or from measuring the flux in the 
UV continuum a t 1500-2800 A in the rest frame of the galaxy 
jKennicutll998l) . The faintness of GRB hosts presents a prob- 
lem for spectroscopy as they require long integration times on 
the largest telescopes. Ground-based photometry in several fil- 
ters presents an alternative possibility for studying the SEDs 
of faint hosts. From such SEDs the UV continuum flux can 
be determined and the galaxy type can be inferred. Previous 
investigations have sho wn that GRB hosts have SEDs simi- 
lar to starburst galaxies Jsokolov et al1l20Q it iGorosabel et alJ 
12003 Jbtlchristensen et al.l2004 . and their SFRs inferred from 
optical methods are moder a te < 10 M^yr"' ( Fruchter et al. 
Il999i iBloom et alJ Il998t bioreovski etail Il998l 12001 
Larger SFRs have been reported based on spectroscopic mea- 
surements. The GRB 000418 host has an un-obscured SFR 
of 55 M Q yr _I derived from the [Oil] line flux ( B loom et alJ 
120031) . All the optical methods for determining the SFRs are 
affected by dust extinction in the hosts. Therefore, the opti- 
cally inferred SFRs represent lower limits to the true SFRs. 
Radio and sub-mm data are much less affected by dust ex- 
tinction, and observations of GRB hosts indicate that the un- 
extincted SFRs can be as much as two orders of magn itude 
larger than those derived from optical estim ators ( Ber ger et"all 
200 ll iBerger et alJl2003l [Tanvir et alJl2004 . However, not all 
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hosts have very large SFRs; som e have SFRs < 200 M H yr 
suggested by radio observations JVreeswiik et all200ll) . 

The SFRs of individual GRB hosts published in the litera- 
ture have been argued to be comparable to th ose of o ther high 
redsh ift galaxies selected by optical methods. iDiorgovski et alJ 
(2001b) found that the [O II] equivalent widths of GRB hosts 
are somewhat larger than that of field galaxies at similar red- 
shifts as GRBs. Likewise. iFruchter et alJ (1 1999b) found that a 
sample of three GRB hosts has bluer colours on average than 
field galaxies in the Hubble Deep Field. 

In this paper we present a statistical analysis of the proper- 
ties of GRB hosts as compared with other high redshift galax- 
ies. We present the 10 GRB hosts selected for the analysis in 
Sect. |2 the derived photometric redshifts in Sect. [3J and the 
SED investigations in Sect.@] We estimate the SFRs of the 10 
GRB hosts by computing the rest frame UV flux in Sect. [5] 
Since the absolute luminosity of the hosts vary by a large fac- 
tor we also analyse the specific SFRs normalised by the host 
luminosities. Comparisons with properties of field galaxies se- 
lected from the Hubble Deep Field are presented in Sect.[6] Our 
results are discussed and summarised in Sect.0 

We assume a flat cosmological model with 
Q„, = 0.3, Qa = 0.7 and H = 65 km s" 1 Mpc" 1 . 
The choice of parameters affects the luminosity distance of the 
galaxies, and thereby the derived SFRs. 

2. GRB host galaxy sample 

Our own multi-colour imaging studies of GRB host galax- 

J I I — i-i 

ies have been presented elsewhere ( Gorosabel et al. 2003a b; 

IChristensen et alj|2004l) . The present work is based on a com- 
pilation of photometry already available in the literature (in- 
cluding our own observational work). 

We imposed a magnitude limit to make sure that the hosts 
entering the sample were bright enough to have fairly accu- 
rate multi-colour photometry in at least 5 optical and near-IR 
bands. This implies that the sample is limited by the available 
multi-colour photometry from the literature. However, with a 
maximum magnitude of R = 25.3 the sample is magnitude- 
limited. These criteria also implied that no host of a GRB oc- 
curring after 2002 is included due to poor multi-colour sam- 
pling. We also required that a redshift for the host galaxy 
or the afterglow be known. Finally, we excluded a few GRB 
hosts which had such complex morphologies that the resulting 
SEDs might be dominated by different sub-components at dif- 
ferent wavelengths (such as the GRB 980613, GRB 011121, 
and GRB 011211 hosts). Having a complex morphology, the 
GRB 980613 host shows colour variations in HST/STIS im- 
ages of Dip — mj p > 0.7 between individual components 

(Hiorthetal. 2002). Similar colour variations was found in 

I I i I 

ground based observations by Diorgovski et al. (2003). Such 
variations in colours make any detailed analysis of the overall 
SED subject to great uncertainty in terms of the derived extinc- 
tion and age. Also the GRB 011121 and GRB 011211 hosts 
have complex surroundings, but not much is present ly known 
about the colours of the hosts themselve s ( Garn avich et alJ 
l2003t Idreineret al.ll2003t I.Takobsson et al.ll2003h . We fhere- 
fore chose to exclude these hosts, but note that these systems 



represent likely mergers which may show significant star for- 
mation. Images of the various hosts are presented elsewher e 
jBloom et alJ2002HDiorgovski et all2003llCastro et alJ2003h . 
To summarize, our selection criteria are: 

- Detection in 5 or more filters 

- R < 25.3 

- Known redshift 

- Not very complex morphology 

These criteria limited our investigations to 10 GRB hosts in 
the redshift range 0.433 < z < 2.037 with a mean and median 
redshifts of z — 0.97 and z = 0.85, somewhat smaller than those 
derived from the 35 GRB redshifts measured to date (z = 1 .43 
and z = 1.10) . The sample is presented in Tabled 

Some of the host magnitudes were obtained a few weeks 
after the burst when the optical afterglow could still contami- 
nate the observed flux. In these cases the host magnitudes are 
derived from fits to the light curves of the afterglows. Since the 
light curve of the afterglow can be described by a power law, 
the total flux is given by /(f) = /o x r a + /host; where the first 
term characterizes the fading afterglow. If the light curve is well 
sampled, the flux of the host /host can be estimated. For exam- 
ple, t he GRB 980703 host magnitudes were derived this way in 
IVreeswiik et alJ (fl999). Data obtained more than one year after 
this particular burst gave magnitudes which are consistent with 
those reported in TablefTIl Hollan d et alJl200ll) . As another ex- 
ample, the expected B band magnitude of the GRB 010921 af- 
terglow w ould be 3 ma gnitudes fainter than the host magnitude 
reported in Pr ice et alJ 112002) at the time of the observations 21 
days after the burst. 

Ideally, the magnitudes of a host should be derived using 
one consistent photometric technique for all filters. For exam- 
ple, in the case of aperture photometry the magnitudes should 
be derived using the same aperture. We can not be sure that 
this is the case for the magnitudes given in Tabled In the cases 
where the hosts are more extended than point sources the effect 
should be negligible as long as the authors have applied a large 
enough aperture for deriving the host magnitudes. 

Furthermore, one should note that we have restricted our- 
selves to analysing only the bright end of the luminosity func- 
tion since it is easier to perform multiband observations of the 
brightest hosts. Only little is known about the nature of the 
fainter host galaxie s JBerger et alJl2002l Ijaunsen et al J 120031 
iHiorth et all2003al) . 

3. Photometric redshifts 

The magnitudes in Table ^ were used fo r comparison with 
theor etical galaxy template spectra from iBruzual & Charlol 
This was done usi ng the program HyperZ 2 developed 
bv lBolzonella et alJ ( l2000t) . 

The templates which are used to fit the GRB host magni- 
tudes consist of elliptical, different types of spiral galaxies, ir- 
regular, and starburst spectra at various ages having different 

1 GRB redshifts can be found at Jochen Greiner's web-pages: 
http : / /www . mpe . mpg . de/ ~ j eg/ grbgen . html 

2 http : / /webast . ast . obs-mip . fr/hyperz/ 
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Host 



GRB 970228 



GRB 970508 



GRB 980703 



GRB 990123 



GRB 990712 



GRB 991208 



Filter (Mag) 


Telescope 


Reference 


B > 26.08 


VLT 


oil j_ l /^An 1 \ 

Sokolov et al. (2001) 


V 


= 25.77 ± 0.2 


nj l/o llo 


Oaldllla CL al. \ ZUUU ) 


R c 


= 25.22 + 0.2 


HST/STIS 




I -- 


= 24.4 ± 0.2 


HST/WFPr? 


Frnrhtpr pt al MQQQa^i 

IT! LLCllLCl CL til. \ lyyya) 


H 


= 23.2+0.3 


HST/NTrMOS? 




K 


= 22.6±0.2 


NTRfVKprk T 


Tharv pt al 1*9009^ 


B 


= 25.89+0.19 


RTA 


OUKUIUV CL al. \ Zajxj 1 ) 


V 


= 25.34±0.22 


BTA 




K 


= 25.06 + 0.17 


BTA 




h 


= 24.11 ±0.25 


BTA 




K 


= 22.7 + 0.2 


Keck T/NTRC 


Charv pt al (900?^ 

11 til y cl til. v z,uuz, } 


B 


= 23.40 ±0.12 


BTA 


Sokolov et al. (200 1 ) 


V 


= 23.04±0.08 + 




Vreeswiik et al. (1999) 


R 


= 22.58±0.06 






I = 


= 21.95±0.25 






J-- 


= 20.87±0.11 






H 


= 20.27±0.19 






K 


= 19.62±0.12 






U 


= 23.6±0.15 




Castro-Tirado et al. (1999) 


B 


= 24.23±0.1 1 




V 


= 24.20±0.15 






R 


= 23.77±0.1 






I -- 


= 23.65±0.15 






K 


= 21.7±0.3 


Keck T/NTRC 


Charv et al (2002 s ) 

\— -lltll V *^ L til. V i-v/v7-_ j 


U 


= 23.12±0.05 


ESO-3.6m 


Christensen et al. (2004) 


B- 


= 23.36±0.09 


DK- 1 .5m/DFOSC 


V-- 


= 22.39±0.03 


DK-1.5m/DFOSC 




R= 


= 21.84±0.02 


UK-1 jm/UrL)SL 




1= 


2 1.41 ±0.03 






.1= 


:20.81±0.17 


NTT7SOFT 




H-- 


= 20.25±0.19 


NTT/SOFT 




Ks 


=20.05±0.1 




T e Flnr'h et al I^OO^U 

Ijt 1 11 CL til. \ jLsVJVJ j 


B 


= 25.19±0.17 1 




castro- uraao et at. (zuui ) 


V 


= 24.55±0.16 




R 


= 24.26±0.15 






I = 


= 23.3±0.2 






K 


= 21.7±0.2 


NIRC/Keck I 


Charv et al. (2002) 



Table 1. Magnitudes in the Vega system in various filters for the selected sample of 10 GRB hosts taken from the literature. The 
magnitudes are not corrected for Galactic extinction. ' indicates that the host magnitudes were derived from power law fits to the 
afterglow light curves. 



star-formation histories. The SFRs decrease with time, such 
that SFR oc exp(-f/r), where the characteristic timescale, t in- 
creases along the Hubble sequence. The starburst template is 
created as an instantaneous burst of star formation (t — > 0), 
which instantly uses up all available gas, while an irregu- 
lar Im template has a constant SFR (r — > oo). We used the 
iMiller & ScaloN 19791) initial mass function (I MF) for calculat- 
ing the templates, as well as a Salneterl (ll955l) IMF with stellar 
masses between 0.1 and 125 M Q . The Miller & Scalo IMF pro- 
duces fewer massive stars compared to a Salpeter IMF and is 
flatter below 1 M . The metallicities of the templates were so- 
lar, Z = 0.02. 

Applying an additional extinction term to the templates, the 
fitting allows an estimate of the type of galaxy, age, and the in- 



trinsic extinction (Ay) for the galaxies. In the fit s we used the 
extinc tion curve estimated for starburst galaxies JCalzetti et all 
2000). We also analysed the SED s using othe r extinction 
curves, e.g. from the Milky Way (fSeatonlll979l) . the Large 
Magell anic Clo ud ( Fitz patrickll986l) . and the Small Magellanic 
Cloud dPrevot et al.lll9841) . We found that the derived ages for 
the dominant population of stars and extinctions did not depend 
on the chosen extinction curve. The reason is that actual values 
of the extinctions are quite small (Ay < 1), thus the various 
extinction curves produce minor differences, as quantified in 
Sect.l4~Tl 

A direct application of HyperZ results in the photomet- 
ric redshifts listed in Table |2] In most cases these are consis- 
tent with the spectroscopic ones. The mean value and stan- 
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Host 



Filter (Mag) 



Telescope 



Reference 



GRB 000210 



GRB 000418 



GRB 000926 



GRB 010921 



U = 23.54 ±0.13 
B = 24.40 + 0.13 

V = 24.22± 0.08 
R = 23.46± 0.10 
/ = 22.49± 0.12 
Z = 22.83 +0.28 
Js = 21.98+ 0.10 
H = 21.51+ 0.23 
Ks = 20.94+ 0.14 

U = 23.54 + 0.3 
B = 24.07+0.05 

V = 23.80+0.06 
R = 23.36+0.05 
/ = 22.79+0.05 
Z = 22.46 ±0.1 
Js = 22.27 ±0.1 
Ks = 21.19 ±0.3 

B = 25.49±0.33 + 

V = 25.08±0.06 
Rc = 24.83±0.07 
Ic = 24.59±0.01 
/a«= 24.1*g 

B = 23.42 ± 0.08 

V = 22.32±0.06 
Rc = 21.93±0.09 
Ic = 21.05±0.08 
J = 20.34 ± 0.02 
H = 19.75 ±0.04 
Ks = 19.07 ± 0.04 



ESO-3.6m/EFOSC2 

ESO-3.6m/EFOSC2 

ESO-3.6m/EFOSC2 

VLT/FORS1 

ESO-3.6m/EFOSC2 

DK-1.5m/DFOSC 

VLT/ISAAC 

NTT/SOFI 

VLT/ISAAC 

ESO-3.6m/EFOSC2 

NOT/ALFOSC 

NOT/ALFOSC 

NOT/ALFOSC 

NOT/ALFOSC 

NOT/ALFOSC 

VLT/ISAAC 

VLT/ISAAC 



Keck I/NIRC 
Keck I/NIRC 
Keck I/NIRC 



Gorosabel et al. (2003a) 



Piro et al. (2002) 
Gorosabel et al. (2003a) 



Gorosabel et al. (2003b) 



Castro et al. (2003) 

J. Fynbo (priv. comm.) 
Price et al. (2002) 

Park et al. (2002) 

Price et al. (2002) 



Table 1. Continued. 



dard deviation of (z P h t - Zspec) is -0.04 and 0.21, respectively. 
Only for the GRB 990123 host the redshift estimate is in- 
accurate (Az > 0.3), but taking into account the uncertainty 
of the estimate the difference is only 1.6<r. A reason for the 
relatively large discrepancy could be that the Balmer jump is 
not well sampled. Indeed, very accurate photometric redshifts 
can be determined if the photometric uncertainties are small 
and the Balmer jump is well sampled as is the case for the 
GRB 000418, GRB 000210, and GRB 990712 hosts. 

At higher redshifts the broad band filters cover a nar- 
rower wavelength range of the rest frame spectrum due to 
the factor (1 + z) accounting for the cosmological expansion. 
Therefore, the accuracy of photometric redshift estimates is ex- 
pected to decrease with increasing redshift. A different mea- 
sure of the accuracy can be obtained by calculating the expres- 
sion Az = (Zphot - z sp ec)/(l + Zspec) for the sample of GRB 
hosts. We find a mean value of 0.015 and a standard devia- 
tion of 0.16 using this expression, showing that it is possible 
to determine photometric reds hifts accuratel y for GRB hosts. 
For comparison. lFernandez-Soto et alJd200ll) derived more ac- 
curate photometric redshifts for galaxies observed through 7 
bands in the Hubble Deep Field, having a standard deviation of 
0.065, which is likely due to smaller photometric uncertainties 
for their galaxies. While the average number of bands of ob- 
servations of the GRB host sample is 6.6, 3 of the hosts have 



observations in 5 filters only. Additionally, the standard broad 
band filters used for ground based photometry are not optimally 
tuned to find photometric redshifts. We find no outliers for the 
estimation of photometric redshifts, indicating that this tech- 
nique is robust for estimating redshifts of GRB hosts. 

4. Spectral energy distributions 

While HyperZ was written for obtaining photometric redshifts 
of galaxies in large surveys, it also serves the purpose of find- 
ing the best matching theoretical galaxy template for a given 
set of broad-band observations. In the remainder of this work 
we shall fix the redshifts of the GRB hosts to the values given 
by the spectroscopic measurements. This was done in order to 
optimise the estimates of other output parameters, as explained 
below. 

The best fit is obtained by minimizing the expression 

i V-l / F host,! — kX FtempJ \ 2 

X =/A TE- \ C 1 ) 

\ o-(F hostJ ) ) 

where the sum is to be taken over all filters, i. Fh os t is the 
flux density of the host in the filter i, o"(Fh st,/) is the asso- 
ciated error, and A: is a normalization constant. F (emp j is the 
flux of the template in the filter i, which is calculated using the 
throughput for the given filter and instrument. HyperZ provides 
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Starburst at z=0.695, age=0.064 Gyr, A v =0.07 




I Host of GRB 970228 



observed wavelength (A) 



Starburst at z=0.706, age=0.091 Gyr, A v =0.37 




Host of GRB 991208 



observed wavelength (A) 



Starburst at z=0.835, age=0.091 Gyr, A v =0.1 7 



rr 




Host of GRB 970508 



5.0-10° 1.0-10^ 1.5-10^ 2.0-10^ 

observed wavelength (A) 




Starburst at z=0.846, age=0.181 Gyr, A v =0.06 



3.0 - 
2.5 - 
2.0 ~- 
1 .5 - 
1.0 '- 
0.5 r ; I 



Host of GRB 000210 



5.0-10° 1.0-10^ 1.5-10 . 2.0-1 0^ 

observed wavelength (A) 



Starburst at z=0.966, age=0.091 Gyr, A v =0.38 




2 - 
t t 



Host of GRB 980703 



observed wavelength (A) 



2.0-10^ 2.5-10^ 



3.0 - Starburst at z=1 .1 18, age=0.064 Gyr, A v =0.14 



2.5 - 
2.0 - 
1 .5 - 



0.5 ~-i 



I Host of GRB 000418 



observed wavelength (A) 



Sa at z=1 .600, age=0.01 2 Gyr, A v =0.96 



Iff 



4-J 




I Host of GRB 990123 



5.0*10° 1 .0-10^ 1 .5-10^ 2.0*10^ 

observed wavelength (A) 



Starburst at z=2.036, age=0.091 Gyr, A v =0.02 



Host of GRB 000926 



observed wavelength (A) 



Starburst at z=0.433, age=0.255 Gyr, A v =0.1 5 




Host of GRB 990712 



5.0-10° 1.0-10^ 1.5-10 . 2.0-10^ 

observed wavelength (A) 



Starburst at z=0.451 , age=0.091 Gyr, A v =0.69 



4 



Host of GRB 010921 



1.0-10^ 1.5-10^ 2.0-1 0~ 

observed wavelength (A) 



Fig. 1. Best fits of the 10 GRB host SEDs fitted to synthetic spectra based on a Salpeter IMF as calculated using HyperZ. 
Redshifts, type of the template, and the extinction are given as inserts in the plots. The solid circles show the available photometry 
for each host (see Tabled after correcting for Galactic extinction. The associated horizontal error bars denote the FWHM of the 
filters. 
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a data set of throughputs for many instruments. For the instru- 
ments where we had no knowledge of the throughput (as for 
example for the Russian 6m-BTA telescope), we simply used 
the throughput of the given filter used for the observation, i.e. 
without convolving with the quantum efficiency of the CCD. 
We tested to see if it had any significance using e.g. the VLT 
throughput curve for the given filter and found the effect to 
be negligible. Since the photometric errors are large in these 
cases, the results for the best fit template did not change. In all 
cases, only templates from one single burst of star formation 
were used. This is a simplification since m ore than one popu- 
lation of stars may be pr esent in the hosts lSokolovetal.12 001: 
IChristensen et alJl2004l) . Hereafter the derived stellar proper- 
ties will refer exclusively to the dominant population. 

Detailed analyses of the SEDs of the host galaxies of 
GRB 000 210, GRB 000 4 18, an d GRB 990712 are pre - 
sent ed bv iGor osabel et al. (2003a), Gorosabel et alJ ( 120 03b). 
and IChristensen et al. (2004), respectively. Similar thorough 
individual SED analyses are beyond the scope of this paper. 

Results from the SED fittings are given in TablefJ] Column 
2 gives the measured spectroscopic redshift, which was held 
fixed while running HyperZ. For comparison the unconstrained 
photometric redshifts are listed in column 3. In column 4 the 
best fit templates are given, in column 5 the age of the tem- 
plate, column 6 gives the extinction, and column 7 lists the re- 
duced x 1 P er degree of freedom (d.o.f.) for the best fits. All 
these values are derived using templates constructed using a 
Miller & Scalo IMF. Column 8 and 9 list the extinction and 
^- 2 /d.o.f. for fits using a Salpeter IMF. The age and galaxy type 
do not change choosing a different IMF, and all the observed 
host SEDs are well fit by starburst templates or young star- 
forming galaxy types. 

Generally, the ages of the dominant population of stars are 
smaller than 0.2 Gyr for all GRB hosts, and the extinctions 
found for the systems are relatively small: < Ay < 2. Plots of 
the best fit templates for each host are shown in Fig.^ 

4.1. Ages and extinctions 

To estimate errors on ages and extinctions we analysed the val- 
ues of x 1 for SED fits of various templates while varying the 
extinction. We derived errors of the extinction and age on the 
basis on which fits gave an increase in the x 2 by A(^ 2 ) = 1. 
The host of GRB 000418 and GRB 980703 are treated sepa- 
rately below. 

Using different extinction laws resulted in differences in A v 
smaller than 0. 1, and furthermore, the estimated ages remained 
constant for any applied extinction law. 

Since GRBs are thought to be produced by the most mas- 
sive stars, ages much larger than 10 Myr of the burst popula- 
tion are inconsistent with this hypothesis. However, we can not 
exclude the possibility that two separate bursts of star forma- 
tion would be able to reproduce the SEDs of the hosts. The 
burst template model uses up all gas in the first burst, while 
a more reliable model would likely have t > 0. We have fit 
all host magnitudes to Im models, which generally produced 
acceptable fits Ot" 2 /d.o.f. < 2). However, this was not the case 



for the GRB 990712, GRB 000210, and GRB 000418 hosts, 
which gave ^- 2 /d.o.f. = 14, 24, and 6, respectively. These 
hosts are the only ones which have extensive multiband pho- 
tometric measurements, and their SEDs are well constrained. 
This implies that a continuous SFR is not in agreement with 
the observed SED. The large errors for the photometry of the 
other hosts make their SEDs less well constrained. 

A special comment is n eeded for the host of GRB 000418. 
In IGorosabel e TaTI (T2003hh . the final reported results for the 
SED fits are consistent with those reported here in Table |2] 
However, the actual best fit having the minimum ^ 2 /d.o.f. 
for a Salpeter IMF and a Calzetti extinction law, is obtained 
with a 0.004 Gyr old starburst templ ate wi t h an ex tinction of 
A v = 1.38 (see Table 3 inlGorosabe l et alJ lE)0 3b) N ). We will 
return to the implication of this uncertainty in Sect. 16.41 

As seen in Table |2] the estimated extinction for the GRB 
980703 host changes by a large factor for the two applied 
IMFs. However, the lower limit on the extinction derived for 
the Miller & Scalo IMF fit is consistent with that derived for 
the Salpeter IMF fit. 

4.2. Ages and metallicities 

There is a well known degeneracy between the age of a stel- 
lar population and the metallicity for a given SED. In order to 
quantify how much this degeneracy affects our results we fitted 
the SEDs to starburst templates, created from a Salpeter IMF, 
with metallicities of 1, 0.4, and 0.2 times solar, respectively, 
using the GALAXEV libra ry of evolutionary s tellar popula- 
tion synthesis models jBruzual & Charloll2003l) . The ages of 
the templates were between 5 Myr and 200 Myr. While keep- 
ing the extinction values fixed to those obtained in Table[2]the 
best fit ages for the different metallicities are listed in Table [3] 
Compared to the ages found by HyperZ in Table |2] there are 
only small differences which are likely due to the differences 
for the input templates. Moreover, there is a general agree- 
ment between the best fit ages for the various metallicities. 
Because of this result we estimate that for the GRB hosts the 
age-metallicity degeneracy produces small systematic errors. 

Alternatively, as shown from the a nalysis of the 
GRB 000418 host in lGorosabel et alJ l2003bl) . the metallicity 
is not strongly constrained in the case where the SED is anal- 
ysed through broad band magnitudes. 

For the median redshift z ~ 1 of the GRB hosts our SED 
analysis covers the rest-frame far-UV to rest-frame J bands. 
Using SED ana lyses of star c l usters in the nearby NGC 3310 
starburst galaxy de Grii s et alJ ( l200 3) find that the ages are well 
reproduced when UV-near-IR pass bands are included, which 
supports our finding of a consistent best fit age. 



5. Star-formation rates 

Since mostly young stars contribute to the UV flux in a galaxy 
there is a relation between the rest frame UV continuum flux 
of a galaxy and the unobscured SFR. One can estimate the 
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Miller & Scalo IMF 



Salpeter IMF 



Host 



^phot 



spectrum 



age (Gyr) 



Ay 



x 2 /d.o.f. 



Ay 



GRB 
GRB 
GRB 
GRB 
GRB 
GRB 
GRB 
GRB 
GRB 
GRB 



970228 
970508 
980703 
990123 
990712 
991208 
000210 
000418 
000926 
010921 



0.695 
0.835 
0.966 
1.600 
0.433 
0.706 
0.846 
1.118 
2.037 
0.451 



0.78*$ 
0.82!™ 

n 4? +0.003 
"■^ -0.01 

0.78*°* 
0.84*$ 
1 00 +002 

1 -""-0.003 

2 09* 54 
37+ 02 

J -0.02 



burst 
burst 
burst 

Sa 
burst 
burst 
burst 
burst 
burst 
burst 



0.064+0.012 
0.091±0.090 
0.006+0.122 
0.012+0.006 

0.255 
0.091+0.090 

0.181 
0.064+0.027 
0.091+0.060 
0.091+0.090 



12* 18 

17* u 

"- 1 '-0.17 

1 60 +005 

-1.27 

90 +017 

"• 7 "-0.20 

15 +004 

"■ 1J -0.05 

35+ 018 

-0.35 

02 +0 06 

"•"^-0.02 

12 +L5 ° 

"• iz, -0.12 

o.oo* 44 

0.68*05 



0.228 
0.146 
0.683 
0.823 
2.815 
0.621 
1.225 
2.448 
0.020 
2.361 



07* 16 

"•" ' -0.07 

17*" 
0.38*^ 

0.96* : ^ 
37+° 16 

"■ J '-0.30 

06* 08 

"■""-0.06 

14 +140 
02* 43 

"■" -0.02 

69* 05 



^/d.o.f. 



0.239 
0.082 
0.831 
0.773 
0.959 
0.613 
1.105 
2.266 
0.025 
2.353 



Table 2. Results from the best fits from HyperZ. Column 2 gives the spectroscopic redshift of the hosts, which was held fixed 
while running HyperZ. Column 3 lists the unconstrained photometric redshifts and the associated 68% confidence levels. The best 
fit template is given in column 4, and the corresponding ages and internal extinctions are given in column 5 and 6, respectively, 
using a Miller & Scalo IMF. Column 8 and 9 list the results of the extinction and^ 2 /d.o.f. for a Salpeter IMF. 



SFR from the re st frame flux at 2800 A using the relation in 
lKennicuttl(ll998l) 

SFR(M yr 1 ) = 1.4 x 10" 28 L V , UV [erg s _1 Hz -1 ]. (2) 

This method is appropriate for obtaining the SFR as long as 
there is a continuous formation of stars. It is a good estimator 
for ages larger than 10 8 years, but shows limitations for ages 
smaller than 10 7 years. Having found that the ages of some of 
the GRB hosts could be smaller than 10 7 years, the relation 
10 is not always applicable. For younger populations the con- 
stant linking the SFR with the luminosity is significantly lower, 
yielding a smaller SFR for a given flux. On the other hand, dust 
is expected to be present in star-forming regions in which case 
the observed flux must be corrected for dust extinction. These 
two effects have an opposite impact on the derived SFRs. 

The method applied to calculate the SFRs from the ob- 
served broad band magnitudes was as follows. First the mag- 
nitudes were corrected for G alactic extinction using the dust 
maps of Schlegel et al. ( 1998). Magnitude offsets calculated by 
HyperZ from standard filters to the AB system were added. 
The AB system is defined as ;«ab = -2.5 log/ v - 48.6 where 
f v is the flux density measured in erg cirr 2 s _1 Hz -1 . Broad 
band magnitudes were converted to flux units (in yuJy) using 
the expression f v = io~°- 4( " ,AB ~ 23 - 9) . The fluxes at the observed 
wavelengths 2800(1 +z) A were estimated assuming power law 
spectra, f v = fo(v/vof for the hosts, between the two filters 
bracketing the observed wavelength at 2800(1 + z) A. 

The calculated SFRs for the 10 GRB hosts are listed in 
Table 0] Knowing the extinction of the hosts from the SED 
analyses, one can correct the SFRs for the effects of extinction. 
Using the extinction curve of ICalzetti et alJ 112000) we derive 
the unextincted SFRs given in column 6. 

The errors of the SFRs in Tableware due to the photometric 
errors used for interpolation which translates into errors of the 
2800 A fluxes. The intrinsic scatter in the calibra tion conv erting 
UV flux into SFR is of the order of 30% JKennicutJl998l) . This 
uncertainty is not included in the quoted SFRs. 



Host 


A v (fixed) 


Z=0.02 


Z=0.008 


Z=0.004 


GRB 970228 


0.07 


100 


50 


50 


GRB 970508 


0.17 


100 


100 


100 


GRB 980703 


0.38 


100 


100 


100 


GRB 990123 


0.96 


5 


5 


5 


GRB 990712 


0.15 


200 


200 


200 


GRB 991208 


0.37 


100 


100 


100 


GRB 000210 


0.06 


200 


200 


200 


GRB 000418 


0.14 


70 


100 


100 


GRB 000926 


0.02 


100 


100 


100 


GRB 010921 


0.69 


100 


100 


100 



Table 3. Ages in Myr derived for fits to various metallicities. In 
addition to the redshifts, also extinction values are held fixed. 



For comparison, SFRs derived from spectroscopic mea- 
surements of the [Oil] line flux are listed in column 7 in 
Table E| Apart from the GRB 000418 and GRB 991208 hosts 
the agreement between the unextincted, UV based and [O II] 
based SFRs is rather good which suggests that extinction does 
not play a major role. Values of the SFR derived either from 
sub-mm, or radio observations are also listed in Table|4] These 
measurements generally show larger values, which may indi- 
cate very obscured components with no (or faint) optical emis- 
sion within the galaxies. 

One immediately sees from Table |4] that the SFRs 
of the hosts are moder at e, in line with t he co n clusion s 
of iFruchter et alJ jl999a|) : iDioreovski etall J2003L l2001al) : 
ICastro-Tirado et alJ feooil) . In Fig. |2] the calculated SFRs are 
plotted as a function of redshift. Since the faint, high-redshift 
hosts do not have multiband observations, and therefore are not 
studied here, the trend for larger SFRs at high redshifts may be 
caused by the selection of intrinsically bright hosts. 

For the GRB 000926 host at z = 2 .037 and the GRB 990 1 23 
host at z — 1.600 we use the observed B and U band mag- 
nitudes, respectively, as a rough measure of their continuum 
flux at 1500 A in the rest frame to get an independent estimate 
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Host 


z 


d L 


UV flux 


SFR 


unextincted SFR 


[0 11 J based SFR (rets.) 


radio/sub-mm SFR 






(cm) 


0"Jy) 


(Mq yr -1 ) 


(Mq yr" 1 ) 


(M G yr~') 


(Mg yr 1 ) 


GRB 970228 


0.695 


1.40X10 28 


0.34+0.16 


0.70±0.32 


0.78±0.36 


0.76(1) 


<335 (8) 


GRB 970508 


0.835 


1.76xl0 28 


0.28±0.15 


0.83+0.45 


1.10±0.60 


1.4 (2) 


<380 (8) 


GRB 980703 


0.966 


2.10xl0 28 


3.20±0.08 


12.7±0.32 


23.8±0.60 


20 (3) 


180±25 (8) 


GRB 990123 


1.600 


3.93xl0 28 


0.55±0.16 


5.71+1.69 


28.0±8.29 






GRB 990712 


0.433 


7.95xl0 27 


1.65±0.04 


1.28±0.03 


1.64+0.04 


2.12±0.60(4) 


<100 (9) 


GRB 991208 


0.706 


1.43xl0 28 


0.40±0.13 


0.83±0.28 


1.54±0.52 


4.8±0.2 (5) 


70±30 (8) 


GRB 000210 


0.846 


1.77X10 28 


0.79±0.07 


2.37±0.20 


2.62±0.28 


3(6) 


90±45 (8) 


GRB 000418 


1.118 


2.52xl0 28 


1.33+0.04 


7.02+0.22 


8.85+0.22 


55 (7) 


330±75 (8) 


GRB 000926 


2.037 


5.28xl0 28 


0.50±0.26 


8.09±4.16 


8.37±4.30 




820+340 (8) 


GRB 010921 


0.451 


8.32xl0 27 


2.15±0.08 


1.81±0.07 


5.67±0.22 







Table 4. SFRs calculated from ground based GRB host observations. Column 2 lists the spectroscopic redshifts and column 3 
the corresponding luminosity distances. Column 4 and 5 list the inferred UV flux and SFRs respectively, and column 6 the SFRs 
corrected for the internal extinction from Table|2] The intrinsic scatter of 30% for the UV to SFR calibration has not been included 
in the reported errors. The last two columns list the SFRs de rived from sp ectroscopic m easuremen t s of th e [O III lines and sub- 
mm/ra dio observations taken fro m fo llowing references: ( l)lBloom et a l] (|2001 ). (2)lBloo m et alJ Jl998l). (3)lDiorgo vski et al 



_ nces: ( 1)l _ . _ . 

( 1998), (4) lHiorthetal](l2000blal) (5) ICastro-Tirado et alJll200ll) . (6) IPiro et alJlhooi . (7) lBloom et alJll2003l) . (8) lBerger et al 



(2003), (9) Vre eswiik et alJ (12001). Note that these SFRs are sometimes derived using a different cosmology than adopted in this 
paper which will cause a small change in the derived SFR. 



15 r 



10- 



LL 

CO 



• f i 



0.0 



0.5 



1.0 1.5 
redshift 



2.0 



2.5 



Fig. 2. Star-formation rates of the 10 GRB hosts as a function 
of redshifts. The SFRs have not been corrected for the effect of 
host galaxy extinction. A trend of larger SFRs for the hosts at 
larger redshifts is very likely a selection effect. 

of the SFR (uncorr ected for extinction). Using the relation in 



of the SFR (uncorn 
iMadau et alldl99Hh 

SFR (M yr" 1 ) 



1.3 x 10~ 28 L V 



[erg s 1 Hz ] 



(3) 



we find SFR = 4.3+2.0 and 7.4+3.0 M yr 1 not corrected for 
internal extinction for the two hosts, respectively. In principle, 
the SFR derived from the continuum at smaller wavelengths 
could be used to constrain the intrinsic extinction, but in prac- 
tice this is difficult given the photometric uncertainties. Indeed, 
one measurement appears to give a smaller SFR based on the 
1500 A calibration compared to the 2800 A calibration, while 
the other is slightly larger. However, in both cases the SFRs are 
consistent within lcr errors. 



5.1. Specific SFRs 

Even though the SFRs appear small, they show variations of 
more than a factor of 10 between the individual hosts. A more 
informative measure of the star-formation activity in the galax- 
ies may be the SFR per unit luminosity. The rest-frame B band 
luminosities are calculated from the best fit (rest frame) spectra 
by convolving with the B band filter transmission. This way, all 
information from the broad band observations is used, and no 
K-correction is involved. Using the cosmological model we can 
estimate the absolute B band magnitude and the corresponding 
luminosity of the galaxy, L, given in Tabled 

It is evident that some of the hosts are less luminous than 
an M* galaxy, i.e., the magnitu de of a galaxy at the break in the 
Schechter luminosity function ( Schec hteJl97(Tl) . A similar con- 
cl usion based on colours o f a sample of GRB hosts was reached 
bv lLe Floc'h et alJ ( 120031) . We adopt M* = -21 which is typi- 
cally inferred for field galaxies. At higher redshifts this magni- 
tude is reported to vary for blue galaxies fcillv et all 1 995). but 
since all galaxies are treated identically here the actual value 
of M* will just introduce a systematic change of the specific 
SFRs. 

The SFRs from Table0]were divided by the quantity L/L* 
in order to calculate specific SFRs presented in Table [5] The 
specific SFRs (not corrected for extinction) vary by a factor of 
~2, as shown in the upper panel of Fig. [3] Choosing a differ- 
ent magnitude for an M* galaxy does not change this result. 
The distribution of specific SFRs has a mean of 9.7 M Q yr 
{L/L*y l and a standard deviation of 2.1. Since the SFR is a 
measure of the flux in the UV rest frame, and the luminosity is 
a measure of the rest frame flux at ~4400 A, the specific SFR is 
simply characterizing the slope of the spectrum for each host. 
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Host 


z 


M B 


LjL* 


specific SFR 










(Me yr (L/L ) ') 


GRB 970228 


0.695 


-18.08 


0.07 


10.3+4.71 


GRB 970508 


0.835 


-18.14 


0.07 


11.6+6.27 


GRB 980703 


0.966 


-21.39 


1.43 


8.87+0.22 


GRB 990123 


1.600 


-21.05 


1.05 


5.45+1.61 


GRB 990712 


0.433 


-18.56 


0.11 


12.1+0.28 


GRB 991208 


0.706 


-18.48 


0.10 


8.45+2.85 


GRB 000210 


0.846 


-19.36 


0.22 


10.7+0.91 


GRB 000418 


1.118 


-20.41 


0.58 


12.1+0.38 


GRB 000926 


2.037 


-20.82 


0.85 


9.55+4.91 


GRB 010921 


0.451 


-19.42 


0.23 


7.83+0.35 



Table 5. Absolute magnitudes, luminosities relative to an L* 
galaxy, and specific SFRs of the GRB hosts. The specific SFRs 
are not corrected for internal extinction. 



6. Comparison with field galaxies 

iLe Floc'h etaP J2003I) found that GRB hosts are rather blue 
compared to other high redshift galaxies, which could be inter- 
preted as GRB hosts having higher SFRs, or being less dusty 
than the average galaxy. We will now investigate whether the 
GRB host SFRs are different from another sample of high red- 
shift galaxies. The SFRs for the field galaxies should be derived 
in the same manner as for the GRB hosts. We therefore need 
a large sample of high redshift galaxies selected from optical 
methods and for which estimates of the redshifts exist. 

6.1. Comparison sample 

Magnitudes in the UBVIJHK bands and photometric redshifts 
of 1067 galaxies in t he Hubble Deep Field North (HDFN) 
have been estimated in iFernandez-Soto et alJ fl999) 3 . Data in 
the UBVI bands were obtained with the WFPC2 using the 
broad band filters F300W, F450W, F606W, and F814W, re- 
spectively. The JHK data were from ground based photometry 
obtained with the 4m telescope at the Kitt Peak Observatory. 
Photometric redshifts are unce rtain within Az=0.1 fo r the 
brightest galaxies with / < 25 dBolzonellaetalJ EoOO). esti- 
mated by comparing the photometric redshifts to the spectro- 
scopic redshifts of more than 100 galaxies in the catalog. Since 
the galaxies have been observed in 7 bands and have smaller 
photometric errors the photometric redshift accuracy is better 
than for the GRB host sample. Additionally, a catalog contain- 
ing 1611 galaxies, with optical data from WFPC2 and near- 
IR data fr om the VLT/ISAAC of the HDF South (HDFS) was 
included dVanzellaetalJl200lh . Compared to the GRB hosts 
these galaxies, which we will collectively refer to as HDF 
galaxies, have a wider span in magnitudes and redshifts. 

The flux densities in the various filters given in the catalogs 
were converted into AB magnitudes which were used as input 
for HyperZ. Conversion factors were calculated using informa- 
tion of the throughputs of the WFPC2 filters for the optical 

3 http://bat.phys.unsw.edu.au/~fsoto/hdfcat.html 



data. Similarly, for the near-IR data the throughputs of the Kitt 
Peak and VLT/ISAAC filters were used, respectively. 

Photometric redshifts, best fitting templates, extinctions, 
and absolute B band magnitudes were estimated for all galaxies 
using the same cosmology as for the GRB hosts. 

In the redshift range corresponding to the redshift distribu- 
tion of the GRB hosts analysed here, 0.4 < z < 2.1, the number 
of galaxies that have photometric redshifts differing by more 
than 0.2 from the spectroscopic measure is 10%. This number 
represents a measure of the overall accuracy of the photometric 
redshift estimations. 

The SFRs of the 1067+161 1 HDF galaxies were calculated 
in the same way as described above for the GRB hosts. In to- 
tal 1140 galaxies were categorized as starburst galaxies, and 
for those the specific SFRs (in M yr~' (L/L*)~ l ) were deter- 
mined. As above, we assumed M* - - 21 and did not correct 
for the effect of extinction. The results are presented in Fig. [3] 
The top panel shows the specific SFRs of 1 140 HDF starburst 
galaxies as a function of their estimated photometric redshifts. 
The specific SFRs (not corrected for extinction) for the GRB 
hosts are shown as large circles. Error bars are not included, 
but for each point the error is -30%, due to the intrinsic scatter 
of the SFR estimator. All specific SFRs of the HDF galaxies 
are in the range 0-20 M Q yr~' (L/L*y l . However relatively 
few galaxies (20%) have specific SFRs exceeding 10 M Q yr~' 
(L/Vy 1 , whereas this is the case for 50% of the GRB hosts. 

To perform a direct comparison with the sample of GRB 
hosts, 851 HDF starburst galaxies with photometric redshifts 
in the range 0.4 < z < 2.1 were selected. The lower panel in 
Fig-E]shows the cumulative distribution of the specific SFRs of 
these galaxies. On average the specific SFRs for the GRB hosts 
are larger than for the field galaxies. Out of the 85 1 HDF galax- 
ies 573 galaxies have specific SFRs above 5 M Q yr _1 (L/L*)~ l , 
which is the lower range of the SFRs of the GRB hosts. Taken 
at face value, this implies that the population of all GRB hosts 
lie among the upper 67% of starburst galaxies, implying that 
GRB hosts have specific SFRs which are larger than for ordi- 
nary field starburst galaxies at similar redshifts. Another ex- 
planation could be that the GRB hosts have less extinction. 
We shall return to a discussion of this effect in Sect. 16.41 and 
Sect.lO 

We performed a statistical test in order to determine 
whether the distribution of specific SFRs vs. redshifts is differ- 
ent for the two samples. We applied a two-di mensional, two- 
sample Kolmogorov-Smirnov test (K-S test) dPeacockl ll983) 
which uses two 2-dimensional samples and checks the proba- 
bility that one sample has the same parent distribution as the 
other. We used the im plementation of the test described in 
iFasano & Franceschinil dl987l). which uses si gnificantly fewer 
computations. According to Peacock ( 1983) the test is valid 
when both the sample sizes are great er than or equal to 10. 
Applic atio ns of the tests, d escribed in IFasano & Franceschinil 
( 1987b and lPeacockl dl983l) . have shown that there is no dif- 
ference for uncorrelated distributions, within statistical uncer- 
tainties, between the two tests. When the probability is > 0.2 
the value of the probability is not accurate, but the hypothe- 
sis that the two distributions are not significantly different is 
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Fig. 4. Age distribution of HDF galaxies and GRB host. The 
solid line shows the ages for all types of HDF galaxies in the 
redshift range 0.4 < z < 2.1, while the dotted line shows the 
distribution of starburst galaxies only. The step-like appearance 
is due to the grid of ages of the available templates. AID K- 
S test give a probability of 48% that the GRB host and HDF 
starburst galaxy distributions are the same. 



0.2 



0.0 — 
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HDF starburst galaxies 
(0.4<z<2.1) 

GRB hosts 



SFR [M yr" 1 (L/L*)" 1 ] 



Fig. 3. Upper panel: Specific SFRs for 1140 starburst galaxies 
in the HDF North and South. The galaxies have photometric 
redshifts estimated by HyperZ and are all classified as star- 
burst galaxies. The SFRs of the 10 GRB hosts are shown as 
large circles. None of the SFRs are corrected for intrinsic ex- 
tinction. Lower panel: cumulative distribution of the specific 
SFR for 851 starburst galaxies having 0.4 < z < 2.1 (solid 
curve) and for 10 GRB hosts galaxies in the same redshift range 
(dotted curve). 67% of the HDF galaxies classified as starburst 
galaxies (573 among 851) have specific SFRs of more than 
5 M yr _1 (L/L*) _1 , which is the smallest specific SFR of all 
the GRB hosts. 



correct JPress et alJll992h . and the derived probabilities can be 
considered as lower limits. 

Our qualitative finding that the distributions of specific 
SFRs vs. redshifts are different for the GRB hosts and HDF 
starburst galaxies is supported by the 2D K-S test, which gives 
a probability of 0.003 that the parent population is the same 
for the two samples. Hereafter, when a two dimensional test is 
performed, the first dimension corresponds to the redshift and 
the second to the SFRs as in the upper panel in Fig. [3] Thus, 
we calculate 2D probabilities for the distributions while show- 
ing the corresponding ID cumulative distribution, e.g. as in the 
lower panel in Fig.|3] 



6.2. Population ages 

In Sect. 0] we found that GRB hosts are young starburst galax- 
ies. We therefore investigated whether they are younger on av- 
erage than field starburst galaxies. For comparison we used the 
ages for the 851 HDF starburst galaxies. In Fig.0]the distribu- 
tion of ages of all HDF galaxies, the GRB hosts and the HDF 
starburst galaxies are shown by the solid, dashed, and dotted 
lines, respectively. 

In addition to the two-dimensional, two sample test, we 
also use a one-dimensional, two sample K-S test. For the ID 
two sample test, the probabilities are reliable for sample sizes 
N = N1N2KN1 + N2) > 4, where Ni and N2 are the number of 
objects in sample 1 and 2, respectively. This criterion is always 
satisfied for the tests performed in this paper. Additionally, the 
probabilities are reliable in contrast with the 2D test. 

Using a one-dimensional K-S test on the distribution of the 
85 1 starburst galaxies ages, we find a probability of 0.48 that 
the two distributions are the same, whereas a 2D K-S test gives 
a probability of 0.18 that the distributions of age vs. redshifts 
are similar for the two populations. Therefore, we have no clear 
indication that GRB hosts are on average younger than field 
starburst galaxies. We furthermore checked whether the red- 
shift distribution of the GRB hosts and the HDF starburst galax- 
ies were the same which is confirmed by a ID K-S probability 
of 70%. 

Comparing the GRB host ages with those of all types of 
HDF galaxies at redshifts 0.4 < z < 2.1 we find a ID K-S 
probability of 0.02 that the distributions are similar. A compar- 
ison of the solid and the dashed lines in Fig.|4]shows that GRB 
hosts are indeed younger than a sample of all field galaxies. 
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6.3. Very young galaxies 

The GRB hosts are classified by HyperZ as very young systems 
with ages smaller than 0.2 Gyr. We therefore also selected a 
comparison sample consisting of young (age < 0.2 Gyr) HDF 
starburst galaxies. As expected, this changes the picture since 
the selected HDF galaxies now have larger average SFRs as in- 
dicated in the upper panel in Fig. EI The GRB host SFRs now 
lie among the upper 83% of the HDF galaxies. Applying the 
2D K-S test gives a probability of 0.01 for the same underlying 
distribution. If one compares the GRB hosts to HDF galaxies 
classified as having ages smaller than 0. 1 Gyr, the two popula- 
tions become rather similar as seen in the lower panel in Fig. [5] 
The 2D K-S test gives a probability of 0.20 for the same parent 
distribution. 

These simple tests show that the specific SFRs of GRB 
hosts are larger on the average than a population of starburst 
galaxies at the same redshifts. The HDF galaxies show a wider 
distribution than the GRB hosts with more galaxies at the high 
and low end of the specific SFR distribution, as indicated by 
the cumulative distribution shown in the lower panel of Fig. [5] 

Our finding that GRB hosts have larger specific star- 
formation rates than field galaxies gives observational support 
to the selection criteria o f potential GRB hosts based on numer- 
ical simulations lCourtv e t al. (2004). Their selected GRB hosts 
do not have large absolute SFRs, but have high star formation 
efficiencies in agreement with our results. 

6.4. Comparisons with extinction corrected SFRs 

The SFRs in Table|5]were corrected for extinction and analysed 
in the same manner as above. Similarly, the SFRs for the HDF 
galaxies were corrected using the Ays estimated by HyperZ. 
The absolute B band magnitudes of both galaxy samples were 
also corrected for extinction, using /?44oo = 5.06 for a Calzetti 
extinction curve (A4400 = RmooEb-v)- Fig.[6]shows the cumu- 
lative distributions of these unextincted SFRs, where the selec- 
tion criteria are the same as in Fig.|5] The upper panel shows a 
comparison of GRB hosts with young starburst (ages less than 
0.2 Gyr) HDF galaxies. The 2D K-S test gives a probability of 
0.04 for the same parent distribution. The lower panel shows 
the younger (age <0. 1 Gyr) HDF galaxies, and the 2D K-S test 
gives a probability of 0. 1 1 for the same parent distribution. 

All the calculated probabilities are given in Table [6] using 
either Miller & Scalo or Salpeter based templates. Generally, 
the differences between the two are small. We can therefore 
conclude that GRB hosts are not drawn at random from the av- 
erage field starburst galaxy population and the GRB hosts are 
most likely similar to HDF starburst galaxies with very young 
(f < 0.1 Gyr) stellar populations. This conclusion is indepen- 
dent on the assumed IMF and extinction correction. 

6.5. Intrinsic extinction 

The probabilities for the extinction corrected SFR distributions 
are generally smaller than for the same uncorrected SFR dis- 
tributions which could indicate that the extinctions found for 
the GRB hosts are different from those of young field galax- 
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Fig. 5. These plots are similar to the one in the lower panel in 
Fig. |5J but here only young starburst galaxies have been in- 
cluded. Upper panel: The solid line represents the specific SFR 
distribution of 689 HDF galaxies with ages less than 0.2 Gyr. 
83% of these HDF galaxies have specific SFRs larger than 5 
M Q yr -1 (L/L*)~ l . Applying the 2D K-S test gives a probabil- 
ity of 0.01 for the sample parent distribution. Lower panel: 327 
HDF galaxies with ages less than 0. 1 Gyr for which 94% have 
specific SFRs larger than 5 M Q yr -1 (L/L*)~ l . The two distri- 
butions are now rather similar, which is confirmed by the K-S 
probability of 0.20 for the same parent distribution. 

ies. In Fig. we show the cumulative distributions of the ex- 
tinction values found for the two samples. Extinction values 
for the field galaxies are from the young (age < 0.1 Gyr) 
HDF galaxies. It seems that the distributions are different in 
the sense that GRB hosts have on average smaller extinctions 
compared to the young field galaxies. However, performing 
one-dimensional K-S tests on the distributions of extinction 
values yields probabilities of 0.29 and 0.08, for the same dis- 
tribution, in the case of the extinctions derived from a Miller 
& Scalo IMF and a Salpeter IMF, respectively. Therefore, the 
small GRB host sample does not allow us to determine whether 
or not GRB hosts have smaller intrinsic extinctions than young 
field starburst galaxies. Likewise, comparing GRB host extinc- 
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Age (Gyr) 



M&SIMF SalpeterlMF I M&S IMF SalpeterlMF 



all 

< 0.2 
<0.1 



0.003 (851) 
0.01 (689) 
0.20 (464) 



0.006 (823) 
0.02 (692) 
0.24 (444) 



0.0001 (111) 
0.0008 (85) 
0.003 (46) 



0.0005 (128) 
0.007 (80) 
0.002 (42) 



extinction corrected 



all 

< 0.2 
<0.1 



0.006 (783) 0.001 (794) 
0.04 (642) 0.004 (676) 
0.11 (439) 0.04(440) 



0.05 (164) 
0.06(131) 
0.008 (84) 



0.001 (218) 
0.006 (164) 
0.004 (95) 



Table 6. 2-dimensional K-S probabilities for the GRB hosts having the same parent distribution of specific SFRs vs. redshifts as 
the starburst HDF galaxies at similar redshifts (0.4 < z < 2.1) with the ages listed in column 1. In brackets are given the number 
of galaxies which are used for comparison. Columns 2 and 3 impose no selection criteria on the brightness of the comparison 
galaxies, while column 4 and 5 include an additional criterion, that the HDF galaxies be as bright or brighter than absolute B 
band magnitude of the GRB hosts. 



tions with those derived for all types of field galaxies does not 
allow us to determine whether they have different extinction 
distributions. 

6.6. Additional tests 

As mentioned in Sect.gjthe SED fit of the host of GRB 000418 
could indicate a somewhat larger extinction (Ay = 1.4). We 
therefore investigated the impact of changing this particular 
extinction value on the derived probabilities for the extinction 
corrected specific SFRs. For a Salpeter IMF, the derived 2D 
K-S probabilities for the same parent distribution are 0.03 and 

0. 04 for 0.2 Gyr and 0. 1 Gyr old HDF starburst, respectively, 
while Table|6]gives 0.004 and 0.04. We therefore conclude that 
the uncertainty in the GRB 000418 host extinction has little 
impact on the results. 

We also investigated the possibility that other galaxy types 
might have similar SFRs as the GRB hosts. Selecting young 
HDF galaxies classified as other spectral types than starbursts, 

1. e. irregular galaxies or spirals with t < 0.1 Gyr, gives a 2D 
K-S probability of 0.002 for the same parent distribution. We 
note that for the theoretical templates the differences between 
the templates at such young ages are small. However, compar- 
ing the specific extinction corrected SFRs with those of ellip- 
ticals, without any age constraint, the probability is also small 
(2xl0~ 5 ) which is not surprising as old elliptical galaxies do 
not have any star formation. 

Because of possible differences between bright and faint 
HDF galaxies, an additional brightness requirement on the 
HDF galaxies was therefore imposed. We selected HDF galax- 
ies with absolute B band magnitudes as bright or brighter than 
GRB hosts (M B < -18.08) and the resulting K-S probabilities 
for the various tests are listed in column 4 and 5 in Table [6] 
The probabilities now appear to be much smaller than without 
the brightness selection (apart from the M&S IMF based ex- 
tinction corrected tests). This is not due to the smaller number 
of comparison galaxies, but arises because the distribution of 
specific SFRs for the GRB hosts is narrower than for the HDF 
galaxies as also seen in Figs. |5] and |6l However, we consider 



that an increase in the GRB host sample is necessary before 
this difference in the distributions can be explained. 

7. Discussion and conclusion 

We have constructed a flux-limited sample (R < 25.3) of all 
GRB host galaxies with known redshifts. The sample con- 
sists of 10 galaxies with broad-band magnitudes in more 
than 5 filters obtained from the literature. The GRBs which 
occurred in these galaxies comprise a collection of a dark 
burst (GRB 000210), a dim one (GRB 000418), a bright one 
(GRB 990712), and a very bright one (GRB 990123). 

Comparing the SEDs of the GRB hosts with template spec- 
tra we find that they are young starburst galaxies with mod- 
erate to low extinctions (Ay < 1). Photometric redshifts are 
found to be accurate, with a standard deviation of +0.21 from 
the spectroscopic ones. Accurate photometric redshifts are ob- 
tained provided there is sufficient optical-IR coverage and the 
magnitudes are accurate to the 10-20% level. Through the anal- 
ysis of the SEDs of GRB hosts we have found that it is impor- 
tant to include near-IR magnitudes when estimating the extinc- 
tion, since the effect of extinction is largest in the UV region, 
and a better sampling of the broad band SEDs of the hosts gives 
a more secure estimate of the extinction. 

By comparing extinctions derived from the SEDs with 
those of galaxies in the HDF North and South, we found that 
the intrinsic extinctions of the GRB hosts are small and not sig- 
nificantly different on average from those of either young field 
starburst galaxies or field galaxies in general. 

The coincidence between small values of Ay from the host 
SED and that of the afterglow suggests that we mostly see 
effects of the global extinction in the afterglow light curves. 
Moreover, small extinction values do not exclude the possibil- 
ity that the GRBs themselves are located in more dusty and 
higher density environments, such as embedded in molecular 
clouds which has been suggested through analyses of X-ray af- 
terglows. This was discussed for the case of the GRB 000210 
host galaxy, where a large Nh was inferred (Pir o et al.l l2002h . 
while the galaxy itself shows a small global extinction (see 
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Fig. 6. Cumulative distribution of the extinction corrected spe- 
cific SFRs of the HDF starburst galaxies with 0.4 < z < 2.1. 
The upper panel consists of HDF galaxies with ages smaller 
than 0.2 Gyr and the lower panel of galaxies with ages smaller 
than 0. 1 Gyr. This again shows that GRB hosts are more likely 
to have similar specific SFRs as very young field starburst 
galaxies. The one dimensional K-S test gives the probabilities 
of 0.04 and 0.11 for the same distribution in the two plots, re- 
spectively. 



Table|2j and lGorosabel et all J2003 a)). However, if the regions 
of star formation where GRBs occur are small and not numer- 
ous, this will not have a large effect on the overall SEDs of the 
hosts. 

For all host galaxies the inferred ages are less than 0.2 Gyr 
while most galaxies have even younger populations, t < 0. 1 
Gyr. A comparison of GRB host ages with those of HDF 
galaxies showed that GRB hosts are not significantly younger 
than starburst field galaxies at similar redshifts, but are clearly 
younger than a sample of all types of field galaxies. 

A good sampling of the redshifted Balmer jump/4000 A 
break gives a well determined age for the dominant population 
of stars in the galaxy. With multiband photometry this jump is 
sampled well for all galaxies in the redshift range involved in 
this study which indicates that the ages of the dominant pop- 



Fig. 7. Cumulative distribution of the extinction values found 
for the GRB hosts (dotted line) and young (< 0.1 Gyr) HDF 
field galaxies (solid line). In both samples a low extinction 
Ay < 2 is preferred by the best fits of the SEDs. The two distri- 
butions appear different from each other, where the GRB hosts 
extinctions are on average smaller than for field galaxies, but 
performing ID K-S tests on the distributions gives inconclu- 
sive probabilities. 



ulation of stars are well constrained. Through fits to templates 
of various metallicities we find that different template metal- 
licities give consistent estimates for the derived best fit ages. 
Specifically, for 2 hosts the ages varied by a factor of < 2, 
while for the 8 remaining hosts, the ages were consistent. 

The SFRs of the hosts were calculated from their rest frame 
2800 A flux and was found to vary by more than an order of 
magnitude from host to host. Specific SFRs, obtained by nor- 
malising the SFRs with respect to the luminosities of the hosts 
are more clustered around the mean value ranging from 5 to 
12 Mq yr~' (L/Uy 1 independently of the redshift. Comparing 
these with specific SFRs of high redshift galaxies in the HDF 
we found that the specific SFRs for GRB hosts lie among the 
upper 66% of the specific SFRs for the field galaxies in the 
same redshift range (0.4 < z < 2.1). We performed several two- 
dimensional K-S tests to quantify the comparisons of GRB 
hosts and subsets of the HDF field galaxy sample. We found 
that GRB hosts most likely have specific SFRs similar to very 
young field galaxies with ages less than 0.1 Gyr. Taking ex- 
tinction effects into account does not change this result. The 
inferred young ages of the dominant stellar populations of the 
GRB hosts are in agreement wit h the idea that GRBs are as- 
sociat ed with core collaps e SNe JWooslevil993tlGalama et all 
fl993lHiorthetalJ2003hl). 

We have found that GRB hosts are not younger than 
field starburst galaxies but have similar specific SFRs as the 
youngest starburst galaxies showing that GRB hosts belong to 
a group of very young, actively star forming galaxies. 

The ages are inferred from the size of the Balmer 
jump/4000 A break, but also from the slope of the spectrum, 
while the specific SFRs measure the flux ratio between 4400A 
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and 2800 A in the rest frame. All SED fits were done with tem- 
plates of solar metallicity, which is likely a simplification. 

Some GRB hosts are found to be Lymana emitters indi- 
cating that th ese galaxies contai n only little dust or have low 
metallicities ( Fvnb o et alJ 120031) . Low dust content and low 
metallicity of the environment is also indicated by spectro- 
scopic observ ations of the optical afterglow of GRB 020124 
(Hiort hetal]l2003al) . Low internal metallicity for GRB hosts 
would imply that we observe bluer colours relative to the HDF 
galaxies which give rise to larger specific SFRs and further- 
more, the SED fits would result in a younger age assuming 
solar metallicity. The conclusions that GRB hosts have simi- 
lar ages as field starburst galaxy and yet appear to have larger 
specific SFRs may therefore be consistent. 

A larger sample of GRB hosts can be constructed by ob- 
taining multiband observations of hosts of bursts which have 
occurred within the past two years. This can be used to analyse 
SEDs and infer SFRs from individual galaxies along the lines 
presented in this paper. Moreover, with future space based mis- 
sions, such as Swift, a sample of uniformly selected GRBs with 
sufficiently brights host galaxies is within reach. This would 
allow a detailed quantitative comparison with the properties of 
specific subsamples of optically selected field galaxies at high 
redshift. 
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